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Abstract 

The early identification of children presenting /UK^'^-mutated neuroblastoma, which are associated with resistance to the 
promising ALK inhibitor crizotinib and a marked poorer prognosis, has become a clinical priority. In comparing the radiology 
of the novel TU-ALK F1174L /Th-MYCN and the well-established Th-MYCN genetically-engineered murine models of 
neuroblastoma using MRI, we have identified a marked /\L/( f " 74i "-driven vascular phenotype. We demonstrate that 
quantitation of the transverse relaxation rate R 2 * (s -1 ) using intrinsic susceptibility-MRI under baseline conditions and 
during hyperoxia, can robustly discriminate this differential vascular phenotype, and identify /WyCTV-driven tumors harboring 
the ALK F " 74L mutation with high specificity and selectivity. Intrinsic susceptibility-MRI could thus potentially provide a non- 
invasive and clinically-exploitable method to help identifying children with MVCW-driven neuroblastoma harboring the 
ALK F " 74L mutation at the time of diagnosis. 
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Introduction 

Neuroblastoma arises in the sympathetic nervous system during 
embryogenesis and is the most common extracranial solid tumor 
in children [1]. For the established subset of patients presenting 
with high-risk neuroblastoma, the current portfolio of therapeutic 
options has limited success, with five-year survival rates rarely 
exceeding 40%. 

The poor clinical outcome and aggressive tumor phenotype of 
high-risk neuroblastoma strongly correlates with amplification of 
the proto-oncogene MTCN and enhanced tumor angiogenesis [2] . 
Recently, mutations in the anaplastic lymphoma kinase (ALK) 
tyrosine kinase gene have been identified in — 8— 1 0 % of primary 
neuroblastoma, leading to constitutive activation of the ALK 
protein. The most common and potent ALK mutation, ALK FII74L , 
is associated preferentially with MTCN amplification, a markedly 
poorer prognosis, and confers resistance to the promising ALK 
inhibitor crizotinib [3-5]. 

With crizotinib in pediatric phase I clinical trials, and other ALK 
inhibitor studies in development, a current challenge is to rapidly 



identify upfront children with high-risk ALK mutated or amplified 
neuroblastoma who may benefit from or become resistant to ALK- 
targeted therapy. 

As most pediatric cancers, including neuroblastoma, originate 
from only a few genetic anomalies during development, they are 
amenable to genetically engineered mouse (GEM) modeling 
approaches. GEM models of neuroblastoma, such as the Th- 
MTCN murine model, which develop spontaneous tumors 
mirroring the major pathophysiological, genetic and radiological 
features of high-risk A/l'CjV-amplified childhood neuroblastoma, 
represent clinically-relevant tools for the study of neuroblastoma 
biology and response to novel therapeutics [6-8]. The develop- 
ment of GEM models co-expressing ALK F1174L and MTCN to the 
neural crest, such as the Th-ALlf" 74L /Th-MTCN model, have 
recently been used to demonstrate how the ALK FII74L mutation 
potentiates the oncogenic activity of MTCN, and that ALK F1174L 
acquired resistance to crizotinib can be overcome through 
inhibition of key cellular modulators of n-mjc [9—11]. 

In this study, we hypothesized that the presence of the 
ALK fJ]74L mutation results in a phenotypic difference in hemo- 
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dynamic vasculature in tumors in the Th-ALE Fn74L /Th-MTCN 
model, and which can be evaluated using intrinsic susceptibility 
magnetic resonance imaging (MRI). The aim of our study was to 
demonstrate that quantitation of the transverse relaxation rate, 
R 2 *, and changes in R 2 * induced by breathing 100% oxygen, 
AR2*oxygen-air, could discriminate between tumors arising in Th- 

ALK F " 7 4L /Th-MTCN and Th-MTCN mice, and that intrinsic 
susceptibility MRI could thus potentially provide a non-invasive 
and clinically-translatable method to help identify children, 
presenting MFC/V-amplified neuroblastoma harboring the 
ALK F1174L mutation. 

Materials and Methods 

Ethics statement 

All procedures involving animals were approved by the Institute 
of Cancer Research Animal Ethics Committee and the UK Home 
Office and carried out according to the United Kingdom National 
Cancer Research Institute guidelines for the welfare of animals in 
cancer research [12]. 

Animal models 

The generation of the Th-ALK FJJ74L /Th-MTCN mice has been 
recently described [11]. Th-ALtf 117 4L /Th-MTCN and Th-MTCJV 
mice were identified by analyzing DNA from mice tails using real- 
time quantitative reverse transcription polymerase chain reaction 
(qRT-PCR, Transnetyx Inc., Cordova, Tennessee), and mice with 
tumors were initially identified by palpation. 

Magnetic Resonance Imaging 

All the 'H MRI studk :s were performed on a 7T Bruker 
horizontal bore micro-imaging system (Bruker Instruments, 
Ettiingen, Germany) using a 3 cm birdcage coil. Anesthesia was 
induced by an intraperitoneal 0.1 ml injection of a combination of 
fentanyl citrate (0.315 mg/ml) plus fluanisone (10 mg/ml) (Hyp- 
norm, Janssen Pharmaceutical, Oxford, UK) and midazolam 
(5 mg/ml) (Roche, Welwyn Garden City, UK) and water (1:1:2). 
A nose piece was positioned for oxygen delivery. 

For all the mice, anatomical T2-weighted coronal and 
transverse images were acquired from twenty contiguous 1 mm- 
thick slices through the mouse abdomen, using a rapid acquisition 
with refocused echoes (RARE) sequence with 4 averages of 128 
phase encoding steps over a 3 x3 cm field of view, two echo times 
(TE) of 36 and 132 ms, a repetition time (TR) of 4.5 s and a 
RARE factor of 8. These images were used to determine tumor 
volumes, and for planning the intrinsic-susceptibility MRI 
measurements, which included optimization of the local field 
homogeneity. The baseline transverse relaxation rate R 2 *, 
sensitive to the concentration of paramagnetic species, principally 
deoxyhemoglobin, was quantified in tumors from Th-ALtf 1174L I 
Th-MTGY(n = 23) and Th-MTCN (n = 21) mice, using a multiple 
gradient echo (MGE) sequence. MGE images were acquired from 
three 1 mm thick transverse slices through each tumor, using 8 
averages of 128 phase encoding steps over a 3 x3 cm field of view, 
and an acquisition time of 3 min 20 s. Images were acquired using 
8 echoes spaced 3 ms apart, an initial echo time of 6 ms, a flip 
angle a =45° and a repetition time of 200 ms. Subsequently, 
100% oxygen (BOC Ltd, Guildford, UK) was delivered at a rate of 
2 L/min to Th-ALjf 1174L /Th-MTCJ\ r (n= 10) and Th-MTCJV 
(n=12) mice. After a 5 minutes equilibrium period, identical 
MGE images were acquired whilst the mouse continued to inhale 
oxygen. 

All the MGE data were fitted voxelwise using in-house software 
(Image View, working under IDL, ITT, Boulder, Colorado, USA) 



with a robust Bayesian approach that provided estimates of R 2 * 
and AR 2 * oxygen _ air ( = R 2 *„ xygen — R 2 * a i r ) [13]. 

Histological Assessment 

Formalin fixed paraffin embedded sections from Th-ALlf 1 1 ' 74L I 
Th-MTCJV and Th-MTCJV mice were stained with haematoxylin 
and eosin, and visualized under light microscopy. The extent of 
functional vasculature in tumors from both models was assessed 
using the perfusion marker Hoechst 33342, and quantified as 
fluorescent area fractions (%), as previously described [8] . 

Statistical analysis 

Statistical analysis was performed using GraphPad Prism 5 
(GraphPad Software Inc., La Jolla, USA). The mean values for 
tumor volume, mean of median values for R 2 * and AR 2 * oxygt , n . air , 
and the mean fluorescent area fractions were used for statistical 
analysis. R 2 * and AR 2 * oxygen . air were assumed to be normally 
distributed, which was confirmed using the D'Agostino-Pearson 
omnibus K 2 normality test [14]. Any significant difference in 
tumor volume, R 2 *, AR 2 * oxygen _ air and the fluorescent area 
fractions between Th-ALlf" 7lt /Th-MTCN and Th-MTCN mice 
were identified using Student's 2-tailed unpaired t-test, with a 5% 
level of significance. 

Results 

ALK F1174L mutation induces a differential radiological 
presentation of neuroblastoma 

Anatomical T 2 -weighted RARE MR images revealed solid 
masses within the retroperitoneum in peri-renal and para-spinal 
abdominal regions of both the Th-ALff 11 74L /Th-MTCJV and Th- 
MTCJV mice, typical of the clinical distribution and radiological 
presentation of human neuroblastoma (Figure 1A) [15,16]. All the 
Th-ALff 1174L /Th-MTCN (n = 23) and Th-MTCN (n = 21) mice 
examined presented with abdominal tumors with wide range of 
volumes (525-3400 mm 3 for tumors in Th-ALlf n74L /Th-MTCN 
mice and 335-2450 mm 3 for tumors in Th-AfTCjVmice). Tumors 
in the Th-MTCN mice appeared heterogeneous, with 94% of 
tumors presenting with areas of hypointense signal, consistent with 
previous observations [8]. Tumors in the Th-ALE Fn74L /Th- 
MTCN mice were generally hyperintense and significantly more 
homogeneous, with only 48% of the tumors having hypointense 
regions (p = 0.006, % 2 -test). On T 2 *-weighted images (Figure IB 
and C, and Figure 2A), tumors in the Th-MTCN mice were 
generally hypointense at longer echo times, whereas the signal 
from tumors in the Th-ALlf n74L /Th-MTCN remained hyperin- 
tense, providing a stark positive contrast with the surrounding 
organs for 74% of the Th-ALlf 11 74L /Th-MTCN mice (compared 
to 14% for the Th-MTCN mice, p<0.001, % 2 -test). 

ALK F1174L mutation results in a slower tumor transverse 
relaxation rate, R 2 *, in neuroblastoma 

Parametric maps revealed a heterogeneous distribution of R 2 * 
in tumors in both models (Figure 2B). Quantitation of the 
transverse relaxation R 2 * revealed significantly slower rates in 
tumors in the Th-ALlf 1174L /Th-MTCN cohort, compared to the 
faster R 2 * rate determined in tumors in the Th-MTCN mice 
(Figure 2D). Baseline R 2 * detected tumors harboring the 
AJjf" 74L mutation with a sensitivity of 90% (95% CI: 66.9- 
98.2) and a specificity of 81% (95% CI: 57.4-93.7). 
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Figure 1. Radiological comparison of the lh-ALK Fn74L nh-MYCN and Th-MYCN mice with abdominal neuroblastoma. A) Anatomical 
transverse T 2 -weighted MR images acquired with a rapid acquisition with refocused echoes (RARE) sequence and B) anatomical transverse T 2 *- 
weighted MR images acquired at increasing gradient echo times as indicated, from representative presenting with abdominal neuroblastoma. C) 
Note the rapidly decaying tumor signal intensity in the Th-MYCN mouse, compared to the more sustained tumor signal observed in the Th-ALK FU74L / 
Th-MYCN mouse. 

doi:1 0.1 371 /journal.pone.0092886.g001 



ALK 



mutation is associated with a stark differential 



BOLD MRI response to hyperoxic challenge in 
neuroblastoma 

Continuous inhalation of 100% oxygen resulted in tumors in the 
Th-MTCN mice demonstrating a very strong, heterogeneous 
increase in T 2 *-weighted image signal intensity (blood oxygen 
level dependent (BOLD) effect) (Figure 2A). In contrast, tumors 
from Th-AIjf n74L /Th-MTCN mice showed a negligible BOLD 
effect on T 2 * weighted image intensity with hyperoxia. As with 
baseline R 2 *, parametric maps revealed a heterogeneous distribu- 
tion of AR 2 * oxy g en - a i r in tumors in both models (Figure 2B and C). 



Tumor regions showing a relatively fast baseline R 2 * remained so 
during hyperoxia, whereas regions of relatively slower baseline R 2 * 
typically showed a more marked reduction in R 2 * with 100% 0 2 . 
The significant difference in R 2 * between the two models was lost 
with 100% oxygen challenge (Figure 2E). As a consequence, 
tumors in Th-ALlf 11741 /Th-MTCN mice demonstrated a signif- 



icantly lower absolute value of AR : 



2 oxygen-air 



than the tumors in 



Th-MTCN mice (Figure 2F). Hyperoxia-induced AR 2 * oxygen _ air 
detected tumors harboring the ALE F1174L mutation with a 
sensitivity of 90% (95% CI: 54.1-99.5) and a specificity of 
94.1% (95% CI: 69.2-99.7). 
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Figure 2. Identification of tumors harboring the ALlf 1174L mutation in MKCTV-driven transgenic mice with intrinsic susceptibility MRI. 

A) Anatomical transverse T 2 *-weighted MR images acquired at TE = 6.8 and 15.6 ms from representative Th-/4/.K F " 74L Arh-MYCW and Th-MYCW mice 
with abdominal neuroblastoma during initial air breathing, and at TE = 15.6 ms after 5 minutes of continuous inhalation of 100% oxygen. B) 
Corresponding parametric tumor transverse relaxation rate R 2 * maps calculated during initial air breathing and after 3 minutes of continuous 
inhalation of 100% oxygen. C) Resulting parametric tumor AR 2 * oxygen . air (R 2 * oxyg en — R2*air) maps. D) Tumor R 2 * during initial air breathing, E) tumor 
R 2 * after 5 minutes of breathing 100% oxygen, and F) tumor AR 2 * oxygen . air (R 2 * 0 xygen — R 2 *air) were determined from Th-zUK^'^VTh-MYCW and Th- 
MYCN mice with abdominal neuroblastoma. Individual data points represent the mean of the median values determined from all three imaging slices 
for each animal, as well as the mean ±1 s.e.m, p, Student's 2-tailed unpaired t-test with a 5% level of significance. 
doi:1 0.1 371 /journal. pone.0092886.g002 
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No correlation between tumor R 2 * and AR 2 * oxygs , n _ air , or 
between either R 2 * or AR 2 * oxyg( , n _ air and tumor volume, was 
determined across the cohorts, indicating that the differences in 
intrinsic susceptibility MRI between tumors in the Th-ALE F1174L / 
Th-MYCN and Th-MYCN mice were independent of tumor size. 

ALK F1174L mutation is associated with reduced functional 
vasculature in neuroblastoma 

Gross examination of the tumors in situ, prior to excision, 
revealed an intense dark red coloration across the whole tumor in 
the Th-MYCN mice, in contrast to the pale appearance of tumors 
in the Th-ALff 1174L /Th-MYCN GEM model (Figure 3A). Histo- 
logical examination with H&E staining revealed the presence of 
large hemorrhagic regions filled with stacked erythrocytes in 100% 
of the tumors from Th-MYCN mice (n= 10), but only in 10% the 
Th-ALlf 11741 /Th-MYCN model (n=10, p<0.0001, x 2 "test) 
(Figure 3B and C). Fluorescence microscopy of Hoechst 33342 
uptake revealed homogeneously and well-vascularized tumors in 
both models (Figure 3D). However, quantitation of the fluorescent 
area fractions revealed that tumors in the Th-MYCN mice had 
significantly higher uptake, consistent with the presence of more 
perfused, functional vasculature (Figure 3E). 

Discussion 

In the current study we demonstrate that baseline R 2 * and 
hyperoxia-induced AR 2 * can discriminate a differential hemody- 
namic tumor vascular phenotype between tumors arising in Th- 
ALff 1174L /Th-MYCN and Th-MYCN models of high-risk, MYCN 
over-expressing neuroblastoma. With a sensitivity of 90% and a 
specificity of 81% for baseline R 2 *, and a sensitivity of 90% and a 
specificity of 94% for 100% oxygen-induced AR 2 *, intrinsic 
susceptibility MRI provides a robust method to discriminate and 
identify Th-MYCN transgenic mice harboring the AIjf n74L 
mutation. 

The rapid and completely noninvasive quantitation of tumor 
R 2 * is suitable for the scanning of young children. Furthermore, 
such pediatric imaging sessions are usually performed under 
general anesthesia, providing an opportunity to transiendy perturb 



the oxygen content inhaled by the patient, and enabling 
quantitation of tumor R 2 * under more hyperoxic conditions (from 
21 to 30-40% 0 2 ). Upon successful translation, intrinsic 
susceptibility MRI could provide a rapid method to identify 
children with ^4LA-driven tumors, enabling the stratification of 
children with this ultra high-risk neuroblastoma at the time of 
diagnosis. 

The baseline transverse relaxation rate R 2 * is sensitive to the 
concentration of paramagnetic deoxyhemoglobin in the vascular 
compartment in tissue (BOLD effect). Compared to most normal 
tissues, tumors exhibit relatively fast native R 2 * values, a 
consequence of the high concentration of deoxygenated erythro- 
cytes within the vascular compartment associated with immature, 
irregular and unstable microcirculation [17]. The fast baseline 
tumor R 2 * measured in the Th-MYCN model is consistent with the 
aggregation of deoxygenated erythrocytes, described as blood lakes 
[18], and which are characteristic of childhood neuroblastoma 
[19]. The significantly slower baseline R 2 * and differential intrinsic 
susceptibility MRI presentation of tumors in the Th-ALK P1174L / 
Th-MYCN mice is consistent with the absence of such blood lakes. 

Inhalation of high oxygen content gases results in the rapid re- 
oxygenation of hemoglobin, with paramagnetic deoxyhemoglobin 
within perfused tumor vessels being replaced by diamagnetic 
oxyhemoglobin, and a reduction in R 2 * [20]. Hyperoxia resulted 
in an overall significant reduction in R 2 * of tumors in the Th- 
MYCN mice, whereas the response in tumors in the Th-ALlF 1174L l 
Th-MYCN mice was negligible, suggesting a clear difference in 
hemodynamic functional vasculature between the two GEM 
models. This was corroborated by the significant overall difference 
in Hoechst 33342 uptake. Interestingly, a spatially different 
hyperoxia AR 2 * response was apparent in tumors in the Th- 
MYCNvaice. Tumor regions exhibiting relatively fast baseline R 2 * 
showed a less pronounced response with 100% 0 2 breathing, also 
consistent with the presence of blood lakes which are typically 
disconnected from the perfused vascular network [21]. In contrast, 
tumor regions with relatively slower baseline R 2 * showed a 
marked response to hyperoxia, indicative of functional vasculature. 

Intriguingly, tumors from the Th-ALR 1174 /Th-MYCN mice 
showed no clear AR 2 * response to hyperoxia, despite Hoechst 
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Figure 3. Pathological comparison of tumors from Th-ALK F " 74L nh-MYCN and Th-MYCN mice with abdominal neuroblastoma. A) 

Gross pathology, B) composite images, and C) high magnification (x200) images from hematoxylin and eosin stained sections. Note the presence of 
large hemorrhagic regions filled with aggregated erythrocytes (*, blood lakes) extravasated from blood vessels (arrowed) in the tumor from the Th- 
MYCN mouse. D) Composite fluorescence images of uptake of the perfusion marker Hoechst 33342 into tumors from Th-ALK F " 74L fTh-MYCN and Th- 
MYCN mice with abdominal neuroblastoma. E) Quantitation of Hoechst 33342 uptake revealed significantly lower functionally perfused vasculature in 
tumors of Th-/4/_/( F " 7 '"/rh-/wyC/V mice (n = 5) compared with tumors in Th-MYCN mice (n = 5). Data are mean ±1 s.e.m, p, Student's 2-tailed unpaired 
t-test with a 5% level of significance. 
doi:1 0.1 371/journal.pone.0092886.g003 
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33342 uptake indicating the presence of perfused blood vessels in 
both GEM models. This suggests a difference in vascular 
architecture that precludes erythrocyte delivery, but not plasma 
perfusion, in the tumor vessels in the Th-ALK F1174L /Th-MTCN 
mice, or in oxygen consumption by the tumor cells [17,22]. The 
extensive hemorrhage and blood lakes present in the Th-MTCN 
model are indicative of vessel wall instability due to rapid 
endothelial cell proliferation and defective pericyte coverage 
[23]. This implies that the tumor vasculature in the Th- 
ALK F1174L /Th-MTCN mice may be more mature. 

Given its relationship to blood oxygen saturation and partial 
pressure of oxygen in and around blood vessels, quantitation of 
baseline tumor R 2 * is also being investigated as an imaging 
biomarker of hypoxia [24]. For example, slow baseline R 2 * has 
been shown to correlate with increased hypoxia, determined by 
pimonidazole staining, in chemically-induced rat mammary 
tumors [25]. Conversely, fast baseline R 2 * is associated with 
hypoxia in prostate cancer [26]. In this context, we have recendy 
shown that neuroblastoma from Th-MTCN mice exhibiting fast 
baseline R 2 *, induced by a hemorrhagic phenotype, are relatively 
oxic, as revealed by negligible pimonidazole staining [8]. R 2 * is 
first of all a marker of impaired hemodynamic function and as 
such shares a causal relationship to hypoxia. However this 
relationship is ambivalent as several differential hemodynamic 
phenotypes can in fact lead to hypoxia, while having an opposite 
effect on R 2 * values. In this study the suggested tumor vascular 
phenotype in the Th-ALK F1174L /Th-MTCN model, which hinders 
the delivery of erythrocytes and causes a significandy slower R2*, 
may result in increased tumour hypoxia in the Th-ALlf 1 174L /Th- 
MTCN model compared with the Th-MTCN model. 

Recent reports have implicated a role of the ALK F1174L mutation 
in tumor angiogenesis in neuroblastoma. Selective targeted 
inhibition of ALK resulted in a significant reduction in vascular 
density in xenografts derived from MTCN non-amplified SH- 
SY5Y cells, which harbor the ALK F1174L mutation, accompanied 
with a decrease in vascular endothelial growth factor (VEGF) and 
matrix metalloproteinases (MMPs) [27]. The expression of MMP- 
9 by stromal cells has previously been shown to regulate the 
vascular architecture in a murine orthotopic MlTiV-amplified 
neuroblastoma xenograft model by promoting pericyte recruit- 
ment [28]. Collectively, these studies suggest that ALK and the 
ALK F1174L mutation contribute to tumor vasculogenesis and 
pericyte recruitment via the regulation of VEGF and MMP-9, 
leading to a dense vascular network with smaller, more stable 
vessels. The differential intrinsic susceptibility MRI vascular 
phenotype observed in tumors in the Th-MTCN and Th- 
ALK FII74L /Th-MTCN models reported herein demonstrates an 
important role for ALE FII74L in angiogenesis in AfTGV-overex- 
pressing neuroblastoma in vivo, strongly suggesting a role in 
vasculogenesis. 
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This study has some limitations. The Th-ALK F,,74L mutation 
leads to the constitutive activation of the ALK protein in 
neuroblastoma [29-32], and the vascular phenotype associated 
with amplified or wild-type (WT) ALK expression is currently 
unknown. Amplification of ALK has also been shown to lead to 
activation of the ALK protein with, however, a 1 7-fold reduced 
kinase activity compared to ALK F1174L mutants [29]. In the 
absence of a Th-ALK"' GEM model, it is difficult to conclude if 
intrinsic susceptibility MRI would be solely able to identify 
ALK: -mutated MYCN amplified neuroblastoma, or more 
generally ^4ZA"-driven MYCN amplified neuroblastoma. 

There is a clear need to more deeply interrogate the role of the 
ALK F1174L mutation on vascular morphogenesis and architecture, 
which may be a major determinant of impaired drug delivery and 
a contributing factor to the poor prognosis of children with 
ALK F1 1 74L -mutated MlXiV-amplified neuroblastoma. The altered 
vascular phenotype may also impact on the response to anti- 
vascular therapies, including cediranib, currently being considered 
in clinical trials for the treatment of high-risk MTGV-amplified 
neuroblastoma [33]. Retrospective analysis of historical patholog- 
ical samples for ALK mutations should provide sufficient statistical 
power to shed a light on the vascular phenotype of neuroblastoma 
associated with ALK amplification, and each of the rare ALK 
mutations, including the lethal ALK F1174L mutation. 

Coupled with our recent identification of quantitation of R 2 * as 
a biomarker of treatment response to cediranib in the Th-MTCN 
GEM model [8], the present study reinforces R 2 * as a biomarker 
of vasculogenesis and its response to therapy in these clinically 
relevant GEM models. Furthermore, it provides a strong rationale 
for the evaluation of intrinsic susceptibility MRI for assessing any 
anti-angiogenic effects resulting from successful targeted inhibition 
of ALK signalling in ALK F1174L mutated neuroblastoma, which 
could ultimately used for the assessment of second generation 
ALK inhibitors [34,35]. To conclude, this study has provided a 
strong rationale for the immediate incorporation of intrinsic 
susceptibility MRI into forthcoming imaging-embedded clinical 
trials of next generation ALK inhibitors in ylZX-mutated and - 
amplified neuroblastoma. 
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